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We show that the ghost degrees of freedom of Einstein gravity with a Weyl term can be eliminated 
by a simple mechanism that invokes local Lorentz symmetry breaking. We demonstrate how the 
mechanism works in a cosmological setting. The presence of the Weyl term forces a redefinition of 
the quantum vacuum state of the tensor perturbations. As a consequence the amplitude of their 
spectrum blows up when the Lorentz-violating scale becomes comparable to the Hubble radius. Such 
a behaviour is in sharp contrast to what happens in standard Weyl gravity where the gravitational 
ghosts smoothly damp out the spectrum of primordial gravitational waves. 
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O . I. INTRODUCTION 

(N 

■ Four-dimensional gravity theories based on action integrals which are non-linear in the curvature invariants have 
rj , been part of the landscape of fundamental physics since Weyl [1] introduced them, soon after Einstein's invention of 

■ General Relativity. 

Such theories possess, apart from the einsteinian ones, new degrees of freedom (dofs) because the action contains 
terms which are non-linear in the second derivatives of the metric. It is commonly believed that, except in f{R) 
theories, see e.g. [2], these new dofs will render even flat spacetime unstable, because they are ghosts, that is, dofs 
O . with kinetic terms with the "wrong" sign, which implies that the phase space of the whole system is no longer finite 
and its energy is unbounded from below. The reason why that is considered dangerous is that any kind of interaction 
$_H . term in the action is expected to yield equations of motion whose solutions are erratic, unstable and eventually diverge 

^ i- 

\ The simplest representative of such pathological theories is Weyl gravity whose action is Einstein-Hilbert's supple- 
^ ■ mented by a Weyl-squared termQ 



S[gab] = ^ jd^X^(R-l Cabcd C^'"") ■ (1.1) 



■ It was analysed by Stelle ^ , who showed the existence of ghosts (non-tachyonic if 7 > 0) but did not exhibit explicitly 
_ ^ , their malignancy since he worked at linear level, where the dofs do not interact. 

■ In [H] ; we analysed linear cosmological perturbations generated during inflation in Weyl gravity, restricting however 
] our attention, as in to the case when the Weyl length scale ^/^ is shorter than the Hubble scale. We found that, 

I ■ despite the fact that no interactions between the dofs were taken into account at this linear level, still the scalar modes 
ILJ ■ diverge in the newtonian gauge, although they remain bounded in the comoving slicing. It therefore seems that the 
malignancy of ghosts shows up already at linear level, when the background is richer than Minkowski spacetime. 
Realising a potential difficulty associated with Weyl gravity already at linear level, we present in the present paper 
^ ] a drastic way out of the potential "horrors" of its ghosts. Following the lines of 0] and Q in their covariant version 
of Hofava gravity [9] we introduce a scalar field that breaks local Lorentz covariance (in defining a preferred time 
direction) and couple it to the Weyl tensor in a way that eliminates its ghosts. However Lorentz- violation will imply a 
modification of the dispersion relation of the remaining, einsteinian, dofs. Specialising to an infiationary background 
we study tensor perturbations and their spectrum (vector and scalar perturbations are not affected by the presence 
of our Weyl term) that we compare and contrast to the behaviour of tensor modes and their spectrum in standard 
Weyl gravity (jl.ip . 



Units: c=l;K = 87rG,7 has dimension length^, k has dimension length/mass, h = k. Mpj . Indices a, 6, • ■ ■ run from to 3; i, j, - ■ ■ 
run from 1 to 3; metric g^i, with inverse , determinant g and signature (— , -f , -(-, +) . R'^tcd = dcV^i,d — dd^'^bc + ■ ■ ■ with V^^c the 
Christoifel symbols; R^d = R'^bad\ R = 9"'' Raf, Gab = Rab -\gabR and Cabcd = Rabcd - 5 (Sac Gbd - Qad Gbc - Qbc Gad + 9bd Gac) - 

■f {gac gbd - 9ad9bc) ■ 
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The organisation of the paper is as follows. In section |TT1 we present our model. In section IIIII we show, in a 
cosmological setting, how it eliminates the Weyl ghosts. Section HVl describes the behaviour of tensor perturbations 
on a de Sitter cosmological background. In section |Vl we quantise these tensor perturbations to obtain the spectrum 
of primordial gravitational waves. Section I VII develops the analysis of [6] of tensor perturbations in standard Weyl 
gravity when ghost dofs are present. Section fVIII discusses and summarises our results. 



II. THE MODEL 

Our gravity model is specified by the action 



where 



S[gab,x] = j^ Jd^xV^ {R + 2jCabcdCefghr'l''^l'''u''u'')+S^[gab,X], (2.1) 



'^"'^ and Jab = gab + UaUb. (2.2) 



Note that the Weyl-term is conformally invariant, see appendix |21 where its relation to the Weyl-square scalar is also 
given. For our purpose the specific form of Sy^[gab, x] for the scalar field x need not be specified, but is requested to be 
such that it yields solutions in which the gradient vector daX is everywhere timelike and future-directed. The vector 
field u"^ then determines a preferred time direction and that necessarily implies that the theory breaks local Lorentz 
covariance. Indeed, the spacetime solution is then preferentially foliated by a family of spacelike hypersurfaces {S^} 
on each of which x takes a constant value. Viewed geometrically, Ua is the future-directed unit normal to T,-^ and jab 
is the induced metric on E,^ . The scalar field x can be called a chronon as in [7| . Note that it can act as an inflaton 
and drive cosmological infiation. 

The equations of motion (coms) for the metric extremise the action with respect to perturbations of the metric and 
read 

Gab - J Bab ^ ^iT^b , (2.3) 

where Gab is the Einstein tensor, Bab is an analogue of the Bach tensor whose expression is given in appendix [Bl and 
T^'' = 2(— g)^^/^ SS^/Sgab is the energy-momentum tensor of x ■ The eom for x , on the other hand, is 



6S^ 7 



K 



VaW^O, (2.4) 



where VaVF" , the variational derivative with respect to x of the Weyl part of the action, is also given in appendix IbI 
As one can see from the expression of Bab in appendix [Bl the Einstein equation (j2.3p contains the derivatives of 
the metric up to fourth order. One may hence worry that the theory should be pathological because higher order 
derivatives theories a priori induce ghosts. 

However, as a closer examination shows, see its expression in appendix |BJ Bab contains only time derivatives up to 
second order if m° is timelike. Since ghosts are induced by the presence of time derivatives higher than the second in 
the eoms, our theory may therefore be free of ghosts. Note that our Weyl action does not give rise either to timelike 
derivatives higher than second in the chronon eom (j2.4l) because the divergence term only contains spacelike derivative 
of the spacelike vector Wa ■ The reason why the field equations are not fourth order in time derivatives is due to 
the presence of the timelike vector , which allows to distinguish time and space derivatives, and thus breaks local 
Lorentz covariance. 



III. GHOST ELIMINATION BY LORENTZ VIOLATION 



Let us first show explicitly and in a specific case that the eoms are indeed second order in time. 
We shall consider the example when the constant x surfaces are intrinsically fiat and take a fiat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) spacetime 

gab dx" dx^ = 0^(77) {-dif + (5„ dx' dx^) (3.1) 



as a background solution, where 77 is conformal time. 
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Note that because of conformal flatness of FLRW spacetimes the Weyl-squared term in the action does not affect 
the background Friedmann equation for the scale factor a(ry) since, then, Bab and Wa vanish, see appendix |B] its 
evolution is therefore the same as in Einstein gravity coupled to a scalar field x = xiv) ■ 

Metric perturbations about a flat FLRW spacetime are decomposed into scalar, vector and tensor parts as 

6gab dx" dx^ = [-2A drf + 2{d.,B + B,) drj dx' + (2C 5.,j + 2d^djE + d^Ej + d^E, + h,j) dx' dx'] , (3.2) 

where diB'^ = diE^ = and where dih^^ = h\ = (all spatial indices being raised with S^^). The following variables 
are gauge-invariant llOi]: 

■f = A+{B-E'y+n{B-E'), <i> = C + n{B-E'), ^, = B,-E[, (3.3) 

where the Hubble parameter is defined by "H = a' / a and where hereafter a prime denotes derivative with respect to 
■q. 

The perturbation to second-order of the action of General Relativity with a scalar field can be expressed in terms of 
gauge- invariant perturbations of scalar, vector, and tensor types [T^. The perturbation of our higher-curvature term 
is found to depend on the gauge-invariant tensor and vector variables, hij and '^i , but not on the scalar variables. It 
reads, see [1| and appendix 

(')^W2N,*z] = j dT^d^X ^g^^^Cabcd^^^Cefghl''^ l"^ 1^' W 

d77d3^(l)Q,fco(^'C^^Ma(^) = l (3.4) 

= i y di^d^'x {dkK^ a'^/i"^" + i Avi/, Avi/'^ , 

where A = 9^9' . 

This is where the role of the chronon is unveiled; it was so coupled to the Weyl tensors that the perturbation of the 
action only contains first order time derivatives, whereas the other possible combinations, Cabcd Cefgh 1°''^ 1^^ 7'^® l'^^ 
and Cabcd Cefgh 7^'^ j'^^ u'^ u'^ , are quadratic in second order time derivatives, see appendix [X] 

It is clear that the vector perturbations are nondynamical, as no time derivatives of appear in p.4p . Moreover 
they are constrained to vanish just as in pure Einstein gravity. 

The total action for the remaining, tensor, perturbations is (see [Tol . [Tl| for the obtention of the Einstein contribu- 
tion) : 

St [h,j] = i- Jd7j d^x [a" {h[j h"^ - dkh.j d'^h'^) + dkKj d'^h'^^ , (3.5) 
and the corresponding eoms are 

1 - j /i-j + 2H - Ah,j = . (3.6) 

We hence see explicitly that, as announced, this action contains only Einstein's gravitons as dofs and that the resulting 
eoms are second order in time derivatives. 

Let us now check that these tensorial perturbations are not ghosts. To do so we decompose them, as usual, into 
two orthogonal polarisations and go to Fourier space: 

h.Av,x)= Y: lj^4(k)hliv)e''-\ (3.7) 



A=l,2 



where the two polarisation tensors e^ (fc), such that e<^^ (fc) e^,* (fc) = Sy , are transverse and traceless, e^^^k) — 
e^^i{k) = 0, and where the reality conditions, 

e^,(fc) = ef/(-fc) and /^^(^) = /,^-*(,y) , (3.8) 

are imposed, a star denoting complex conjugation. Then the action p.Sp reads, in Fourier space 

STiH}]-^-. E fdvd'k[ia^+A^k')\h'l\'-aH'm'], (3.9) 
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where k = \k\ . The graviton can thus be viewed as a coUection of non- interacting scalar fields as in General Relativity. 
The important point to note here is that 7 must be positive, otherwise the graviton modes would be tachyonic ghosts 
on Minkowski spacetime when 0(77) — 1. Therefore (a^ +47^^) is always positive, so that the kinetic terms \h'^\'^ in 
p.9p are always positive, and the graviton never becomes a ghost on a FLRW background. 

IV. EVOLUTION OF TENSOR PERTURBATIONS IN INFLATIONARY COSMOLOGY 

What is to be checked now is whether or not the Weyl term modifies drastically the time evolution of the modes. 
The eom in Fourier space deduced from p.6p and p.9[) is: 



x';: + nlxi: = with xj:iv) ^ Vo^+^ h{r,) , (4.1) 
where from now on we omit the index A and where the pulsation ilk{v) is given b}Q 



a2 + 47fc2 (a2+47/s2)2_ 

We shall first discuss the time evolution of the modes on a de Sitter background. Setting 

1 



(4.2) 



and z = —kri, (4-3) 
—H 77 

H being a constant, the eom reduces to 

d^Xfc- , - 2 + 4^2^2 (^2 _ 3) fcVl + 4e2z2 

-d^+ z2 (1+4,2,.). X,-0 w:th x,= h, (4.4) 

where we have introduced the dimensionless parameter 

e^VlH. (4.5) 

Equation (|4.4p can be solved exactly in terms of hypergeometric functions, see next section. It is however enlightening 
to study first the qualitative behaviour of the modes. 

As in General Relativity, a mode with wave number k becomes larger than the horizon at late times, when z <^ 1, 
that is when its physical wavelength a/k — 1/{H z) becomes much larger than the Hubble radius . As for the 
Lorentz violating regime it holds at early times when ez ^ 1, that is when the physical wavelength a/k oi the mode 
k is still much shorter than the length scale ^/j set by the Weyl term. 

Let us first consider the case e <C 1, when the length scale ^/7 on which the Lorentz-violating Weyl term operates 
is much shorter than the Hubble radius H^^ . (This is the case if the Weyl correction is viewed as a low energy limit 
of some Planck scale quantum theory of gravity, inflation happening at, say, the grand unified theory (GUT) scale.) 

The modes Xk ' which solve (|4.4L will then go through three different regimes: the early, Lorentz violating regime, 
when e z 3> 1; an intermediate regime, 1 3> ez 3> e when they have left the Lorentz violating period but are still 
sub-horizon; and the late stage, z ^ 1 when they have exited the Hubble scale. 

In the late and intermediary regimes, when ez -C 1, which are Lorentz symmetric, equation (j4.4p is the same as in 
General Relativity: 

^f+(l-|)«=.0. (4.6) 

whose two independent solutions are well-known: after horizon crossing, that is for z <C 1, we have cx 1/z and 

X^^^ cx z2 , so that the growing tensor perturbation behaves as hj: cx zx^^'' — > const., i.e., "freezes out"; and on 
sub-horizon scales z 3> 1 (but z <^ 1/e) the modes xt oscillate as sinz and cosz . 



The analysis of the modes can also be performed in terms of cosmic time t such that dt = a{ri) drj , in which case the eom reads 
fk + '^'^fk = with fk = a^/^ y'l -I- 47 (k/a)'^ hk , where a dot denotes derivative with respect to cosmic time, where H = a/a and 
where 



4^ ' l + 47(fe/a)2 [l + 47(fc/a)2]: 



5 




FIG. 1. The growing and decaying mode functions x[f''(^) s^nd x^^\^) on a de Sitter background, with z — —kr) = k/(aH) , for 
e = ^ H <^1, that is, when the energy scale k Mpi/ of the Weyl-term is much higher than the inflationary scale k Mpi H . 
In the above figures, where e = 0.01 and time increases from right to left, the Lorentz violating period ends at around z — 100 
and Hubble radius crossing occurs at around z — 1. 



In the, early, Lorentz-violating regime, e z ^ 1 , oc /i^: , and the eom (j4.4p reduces to 

^^Xfc _|_ Xfc ^ Q 



(4.7) 



whose solutions oscillate as xj: c« 2^/^=*="^/^ , with v = -^/l/e^ — 1 real since e < 1. See figure[TJ 

Let us now consider the case e ^ 1, when the length scale on which the Lorentz-violating Weyl term operates is 
much longer than the Hubble radius . (This would be the case if the Weyl correction was emerging from some 
effective theory operating below the inflationary, GUT, energy scale.) 

At late times, e z ^ 1 ^ e the modes behave as before and as in General Relativity: D,'^ ~ — so that x^?"* oc 1/z 
and x^'^^ oc . In the intermediary stage on the other hand, when 1 ^ e 2 <C e , they are already in the Lorentz- 
violating stage and since, then, fll ~ — 3/(4e^z'^) , they behave as Xk °^ ze'^^^'-'^'^^^ and do not oscillate. But the 
really important difference with the previous case occurs at very early times when e z ^ e ^ 1. We still have then that 
rtj: ~ l/(4e^ z^) but the modes, instead of oscillating, behave as Xk ^ z^l'^'^^l'^ with v = yj\ — 1 / real. Because they 
never oscillate we can qualify these modes as "rampant" , that is, as "flourishing and spreading unchecked" (Oxford 
dictionary). 

We have up to now approximated inflation by a de Sitter stage. In the case of slow-roll inflation, the Hubble 
parameter H is no longer constant but slowly decreases with cosmic time, i.e. with decreasing z . The previous results 
then hold when the parameter e — slowly decreases. 

Thus, if the energy scale at which inflation starts is lower than the energy scale set by the Weyl term, in other 
words if the Lorentz-violating length scale is always shorter than the inflationary Hubble radius (e = ^ H < 1), then 
all modes behave as in figure [T] If now inflation starts when the Hubble radius is shorter than the Lorentz-violating 
scale, then the modes whose wavelengths are big enough to exit the Hubble scale when e is still bigger than 1 will 
never oscillate. On the other hand shorter wavelengths modes, which will still be within the Hubble radius when e 
goes through the critical value 1 will start in a non-oscillatory way at early times and then will behave qualitatively 
as in figure [TJ once e becomes smaller than 1. 

Let us summarise sections IIIII and IIVI we confirmed that our Lorentz-violating Weyl action (|2.1[) , when considered 
in a cosmological setting, indeed yields second order equations of motion for the tensor perturbations (which are 
the only ones to be modified by the presence of the Weyl term), see p.6p : we confirmed also that, for 7 > 0, the 
perturbations never become ghost-like despite the modification of the coefficient of the kinetic term of the graviton due 
to the Lorentz-violating Weyl term, see p.9p : finally we studied the perturbation modes, saw that they differ markedly 
from standard Einstein-de Sitter gravitons in the distant past, and distinguished two different early time behaviours 
according to whether the Lorentz-violating scale is shorter or longer than the Hubble scale. These modifications 
of the early time behaviour of the modes will induce a significant change in the spectrum of quantized primordial 
gravitational waves, as we shall now see. 
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V. QUANTISATION OF TENSOR PERTURBATIONS AND PRIMORDIAL GRAVITATIONAL WAVE 

SPECTRUM 



Let us for clarity write again the action (j3.5p , for the tensor perturbations on a FLRW background with scale factor 
a{r]): 

The momentum conjugate to hij is 

tt'J' = — (a^ h"'J - 47 A/i'*^) . (5.1) 
4k 

Canonical quantisation turns hij and tt*^ into operators satisfying the standard commutation relation 

[hij{ri,xi), 7r'-'(77,X2)] = 2ih5{xi - X2) , (5.2) 

all other commutators being zero (the factor 2 comes from the fact that hij is the sum of two independent dofs). 
Expanding hij in Fourier modes as in (j3.7p . we further decompose h^ as 

h\{i^)^alhk{v) + fitihk{^). (5.3) 

where and are arbitrary operators (we introduce - for later convenience) and where h^ and hk (77) are two 

independent solutions, depending on k = \k\ only, of the second order eoms in Fourier space deduced from p.6|) . that 
is: 

1 + /ifc + 2n h'k + k^hk^O, (^1 + h'l + 2H h'k + k^hk = 0- (5.4) 

It follows from the eoms (|5.4p that the Wronskian of hk{ri) and hkii]) j 

W{hk,h) = hk [a" + 47 fc2) /i;^ - hk [a" + 47 fc2) h'k , (5.5) 
is a constant ^ if /ifc (77) and /ifc {rf) are independent solutions of (j5.4p . Hence 

/i..(^,2')-E /(^^ 4(fc)a^;ife(r7)e"=- + 4*(fc)5^-/i,(r;)e-''=-~ , (5.6) 



A=l,2 



where we used the reality condition on the polarisation tensors p.Sp . (Note that we have not yet implemented the 
hermiticity condition on h'^{ri) .) 

Plugging the expansion (|5.6I) into (|5.ip . we find that the commutation relations (|5.2p are equivalent to 

[alralj^Sih-h) (5.7) 

(all other commutators being zero) if the Wronskian of the two independent mode functions hk and hk is normalised 
to 

W{hk,hk) = hk {a^ + A'jk'^)h'k - hk (a"^ + k"^) h'k = Anhi. (5.8) 

In Minkowski space we have a — I and the eoms (|5.4p for the two independent solutions hk (??) and hk {rf) reduce to 

^ - fc2 

hk + i^Mk hk = 0, hk+ WMfc hk=0 with a;^^^ = _^ ^ , (5.9) 

where Wma; positive and where, for clarity, a dot denotes derivation with respect to cosmic time t such that 
dt — a(r]) drj . 

The "positive frequency" modes, which will define the vacuum |0) such that |0) = and at |0) = 0, are chosen 
to be 

/ifc = rjfc e-' /ifc =nfc e+' (5.10) 
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where the coefScients and must be such that the hermiticity and Wronskian conditions, see (I3.8P and (|5.8p . are 
satisfied, that is, such that 



a^Uk^a^'Uf., nknk = p , (S-H) 



which impose (up to irrelevant constants) 



a^, = af, n, = nl^^^^^. (5.12) 



Therefore, the choice for the modes is 



(5.13) 

k 

Note for further reference that in the short wavelength limit, j k'^ — oo , the pulsation becomes independent of k : 
WMfe — >■ (47)^"^^^ , so that we have 

/..^^e--/(^v^). (5.14) 

Next we take accelerated cosmological expansion into account, where another scale H , the Hubble parameter, 
comes into play. 

In the de Sitter case when a — l/{—Hri) with H = const, and rj conformal time, the eoms for the modes as given 
in (j5.4l) read 

9, (Phh 2 dhi. hi. , (Phh 2 dhi. hi. 

where 

e = y/j H and y=—€kr]. (5.16) 
As for the Wronskian normalisation condition (15.81) it becomes 



dhk J dhk . AnhH^ 

Now, two independent solutions of (|5.15p are 

L ,^ I ^ f -l-iv -l + \v I , 2\ , /X 32 3 ^ /5 + i:/ 5-iiy 5 .\ 

^(g)(y) = 2^(^^'^^'-2'^^^ ) ' '^(d)(y) = y2/^(^^'^^'2'~^^ J ^^-^^^ 



with 



|iA/l-l/e2 if e>l ^ ' 

Note that the hypergeometric functions introduced in (15.181) are real, whether v is real or imaginary, and that 

/1(g)(0) = i, /1(d)(0) = 0. (5.20) 

In order now to determine which combinations of ft.(g) and /i(d) we must choose as our independent modes hk and 
hk we have to study the early-time, large k limit. 
When y — ^ 00 we have 

h^^^iy) ^ C(-^^ (2y)V2+i^/2 + (22;)V2-i^/2 ^ ^.^^^(y) ^ (2y)V2+i./2 + ^ (2^)1/2-1^/2 (5 
with, when is real, 

_^¥r(-W2^ ^ ^_^¥IVW2^ 
^(g)-^(g) - r2(-i/4-ij./4) ' V)-^(d) - r2(5/4-iz./4) ■ ^ ^ ^ 

(When h' is imaginary then the coefficients are real but unequal.) 
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In cosmic time t , such that dt = adrj , we have a 



and H y = kee , so that, for v real (i.e. for < 7 H < 1) 



with 



2^7 



which, when ^ <C 1, identifies with the Minkowski positive high frequency modes chosen in (|5.10p : 



y 



l/2+ii//2 ^ g-it/(2V7) 



(5.23) 



(5.24) 



We recover here what we have discussed in the previous section, to wit that the early time behaviour of the modes 
differs drastically depending on the value of j : they oscillate if i' is real (7-ff^ < 1), and are "rampant" if 1/ is 
imaginary (7^?^ > !)■ It is clear that quantisation will make (easy) sense only when modes can qualify as "positive 
frequency" , which implies that they must oscillate. We shall therefore suppose from now on that 



< 7^2 < 1 



0<e= y/jH < 1 



1 real . 



(5.25) 



In keeping to what we chose when the background is Minkowski spacetime, we shall hence impose that the two 
independent solutions hk and hk behave, for large y , as 



hk{y) ^ nu (2y)i/2+i-/2 , -hk{y) -> (2y)i/2--/2 . 



(5.26) 



As for the values of the coefficients nk and rik they follow from the large y limit of Wronskian normalisation condi- 
tion (|5.17p as well as the hermiticity condition, see p.8|) . which imposes, up to irrelevant constant that 



Uk = Uf, 



K h 



2v fc3 e3 



H. 



(5.27) 



Thus, all in all, the hermitian operator hij reads 



A=l,2 



(27r)3/2 



4(fc)a^/ife(,7)e"=- + h.c 



(5.28) 



where the modes are: 



8fc3 £3 



(d) '^(g) "C(g)/i(d)), 



(5.29) 



where the functions /i(g) and ft-(d) a-nd the coefficients C(g) and C(d) are defined in (j5.18p and (j5.22p (and where we used 



-(d) '-(g) '-(d) ^(g) 



= 2i/v). Finally the "Bunch-Davies" vacuum |0) is defined as usual by 
«^|0)=0. 



The power spectrum V{k] rj) of the gravitational waves htj is now defined as 



\h,,iri,x,)h'^{r^,X2)\0) 



d^k 



Ank^ 



and is given by, using the commutation rule (j5.7l) 



Vik;r]) 



fc3 



(5.30) 



(5.31) 



(5.32) 



After horizon crossing, that is in the limit 77 — >■ 0, we have hk oc fc^3/2 ^ ggg (j5,29l) . and the power spectrum is scale 
invariant just like the power spectrum of ordinary gravitational waves in Einstein theory on a de Sitter background. 
However its amplitude depends on 1/ = Vl/e^ - 1 and is, using the late time limit of ft-(g) and ft.(d) , see (|5.20p 



2kH^_ _ _ cosh(7rr//2) coth(7rzy/2)r2(-l/4 + ii//4)r2(-l/4-ii//4) 

I \rh] Tj — > U j — , WLier6 — 



1287r2 e3 



(5.33) 



where we used the relations r{-iiy/2)T{iv/2) = 27r/[i/ sinh(7r j//2)] and r(5/4 - i i//4) r(5/4 -I- i z//4) r(-l/4 - 
it//4)r(-l/4 + ij//4) = 27rVcosh2(7r://2), see figure [3 
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FIG. 2. The modification factor of the power spectrum of the tensor perturbations on a de Sitter background as a function of 



This is the main result of the paper: In our ghost-free but Lorentz violating Weyl theory of gravity the spectrum 
of gravitational waves, when evaluated on a de Sitter background, is scale-invariant, as in General Relativity, but its 
amplitude in the late time limit varies with e = ^ H , where is the length scale on which the Weyl term operates 
and H^^ is the Hubble radius. For small e, S w 1: the spectrum of primordial gravitational waves is almost the 
same as in General Relativity because the Lorentz symmetric regime before horizon crossing lasts long (it spans 
a z- interval much bigger than 1); as e increases, that is as the Lorentz- violating length scale approaches the Hubble 
radius, its amplitude decreases up to a factor 65 %, before eventually blowing up as e grows further and approaches the 
critical value 1. For e > 1 the spectrum, as we saw, cannot be defined as the vacuum expectation value of quantized 
modes since they do not oscillate at early time. Another criterion must then be chosen to define it in this regime. 

In any case, the spectrum of primordial gravitational waves is greatly modified, if the Lorentz- violating scale happens 
to be of the same order of magnitude or larger than the Hubble radius. But, again, it must be stressed that only 
when the Hubble radius is considerably larger than the Weyl-correction characteristic length scale (e ^ 1) do the high 
frequency modes become those of flat spacetime and it is only then that there is a natural definition of the vacuum 
state. 

In slow-roll inflation now, the Hubble parameter H is no longer constant but becomes a slowly decreasing function 
of cosmic time. The previous results therefore still hold but H (and hence e = and v = -^/l/e^ — 1) in the 

expression (j5.33p of the spectrum may be evaluated around horizon crossing when H{t) = k/a{t) . This allows to 
express t , and therefore H , e and in terms of k . Since e(fc) is a decreasing function of k , the qualitative behaviour 
of gravitational wave spectrum V{k; i|ij=/c/a) is then read off from the right to left in figure [2] when k increases. What 
happens when e{k) crosses the critical value 1 during inflation remains to be investigated. 



VI. PRIMORDIAL GRAVITATIONAL WAVES IN EINSTEIN PLUS PURE WEYL SQUARE GRAVITY 



In order to compare and contrast the ghost-free but Lorentz-violating Weyl theory of gravity that we have studied 
in this paper with ordinary "ghastly"0 Weyl theory, we summarise and develop here the results of [6] to obtain the 
exact spectrum of gravitational waves on a de Sitter background in Weyl gravity. 

When expanded to second order in the tensorial perturbations around a FLRW background, see p.2p for definitions, 
the action of Weyl gravity (jl.ip becomes (see see (5j] or appendix E]) 

S[h,,] jdv dPx [a" {h[j h''^ - dkh,, d'^h^') - 7 {h'^^ h"'' - 2dkK^ d'^h'^' + dkiK, d'^'h'^)] , (6.1) 

which must be compared to (|3.5p . (Recall, see [1] , that 7 must be positive to avoid the ghost dofs to be tachyonic on 
Minkowski spacetime.) Proceeding as in 16[ "a la" Ostrogradsky we introduce a new variable Qij = as well as a 



^ i.e. "causing great horror or fear" (Oxford dictionary) 
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Lagrange multiplier A*-' (both of them transverse traceless) and consider the equivalent action 



S[hij , Qij , A ■ 



The conjugate momenta are 



8k J 



— (a^ h"^ -2jAh''^ - X'^) . 
4k 



4k 



The eoms obtained by extremisation of the action with respect to A*-' and Qij are 
and allow to express the momenta in terms of hij alone as 



1 



Quantization is implemented by imposing the commutation relations 



4k 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



\h^j(ri,xi), ttJ;' (77,^2)] = 2ift(5(xi - xi) , [Qij(r], fi), 7r^(?7, ^2)] 2i?i(5(fi - X2) , 



(6.6) 



all other commutators being zero and the factor 2 coming from the fact that both hij and Qij represent two dofs. 
We now expand hij in Fourier modes 



hij{r],x) 



(27r)3/2 



^ e^^. (k) al hi'^ (r;) + ^ ef^ (fc) 6^- /.f (r;) e^-^V h.c 



A=1.2 



A=1.2 



(6.7) 



Plugging this expansion into the definitions of the momenta (j6.5p . we find that the commutation relations (j6.6p are 
equivalent to imposing (omitting the index A) 



[a,-^ , dl ] = <5(fci - fc2 ) , , St ] = -S{h - 



(all other commutators vanishing) if the mode functions satisfy the following Wronskian conditions 

[{a' + 27 fc^) h'^'^*+jhf'^*] - [{a' + 27 fc^) hf^* +lhf'^*] - c.c. = 4Ki?., 



l'(2) , "(2)* 



Anih 
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(6.8) 



(6.9) 



The constancy of the Wronskians is ensured if the mode functions and and their complex conjugates are four 
independent solutions of the eom for hij which reads, in Fourier space 

(a^ h'J + aH^hk+-/ (hi" + 2k^ h'l + k^ hk) = . (6.10) 

Let us now specialise to a de Sitter background where a = l/(—Hri) , H being a constant and 77 conformal time. 
Then, as shown in @ the eom (|6.10p factorises neatly: if /i^^-* = z ii^j^^ and h^j^^ = z /x^j^-* are taken to solve, respectively, 
(with z = —krj) 



dz2 



3.11) 



7i?2 z2 



then, as can be shown explicitly, they satisfy the original eom (|6.10p and, also, the Wronskian conditions 



if 
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The first equation (|6.1ip can be seen as the eom for the usual Einstein graviton and the second for the Weyl ghost dofs. 

What remains to be done is to choose "positive frequency" modes which define a "Bunch-Davies" vacuum state |0) 
such that Gk |0) — bk |0) — 0. This is easy in this case since there is no violation of Lorentz covariance, so that the 
eoms (|6.1ip both reduce to the standard Minkowski eoms for massless fields tor z ^ 1 and z ^ 1/(7 i?^) , that is in 
the remote past when the physical wavelengths a/k of the modes are much shorter than both the Hubble radius scale 
H^^ and the length scale ^7 set by the Weyl term. We therefore impose that for z — > 00 , 

, ^ ze''. (6.13) 

The solutions of the eoms (j6.11l) having this early time/short wavelength behaviour are, up to a phase 




(1) . /ZK H 



hl'>{z)^d—^^^={l + iz)e'\ 



(6.14) 



with D EE (1/2) ^1 -4/(7i?2). 

The power spectrum V{k; z) of the gravitational waves hij is again defined as 



and is given by 



P(fc;.) = ^ (|41)(.)P_|42)(,)|2) 

2kH^ 1 



7r2 l + 27i72 V 2 



l + ,2_ " 3|^i../2^(l)(^)| 



(6.15) 



(6.16) 



Whatever the value of 7ij2 , that is whatever the sign of 1^2 — 1 [1 _ 4/(7ij2)] ^ the last two terms do not contribute 
to the power spectrum when z 0, that is, at late time when the modes have exited the Hubble radius, and we have 
that 

, 2Kij2_ _ 1 

'P[k;z^O) = 5— ^w, where = , , ^ ■ (6.17) 

TT-^ 1 + 27 ii^ 

(This generalises @, where the spectrum was computed for small values of the parameter 7_ff2 only.) 

This de Sitter spectrum of primordial gravitational waves at late times, obtained in Einstein plus pure Weyl-square 
gravity based on the action (|l.ip , and that obtained in (|5.33p and figure [5] in the ghost-free but Lorentz- violating 
Weyl gravity theory based on the action (|2.ip . are very different: contrarily to S, Sw never blows up and, as 7ij2 
tends to infinity, goes to zero as {'f H^)~^ . 



VII. CONCLUSION 



We have proposed and investigated a mechanism to eliminate the ghost degrees of freedom of higher-curvature 
gravity theories. The mechanism invokes the gradient of a scalar field which is timelike and therefore implies a 
breakdown of Lorentz covariance. Although we only discussed here the particular example of Weyl gravity, the 
mechanism appears generic enough to eliminate any higher-derivati ve g hosts. We expect that a canonical way of 
analysing generic higher-curvature gravity, developed by the authors 12], will be useful in studies in this direction. 

We investigated quantisation of the inflationary tensor perturbations in our ghost-free Weyl gravity model and 
defined the "positive frequency" modes as those reducing to flat spacetime, Lorentz-violating, positive frequency 
modes. Such a modiflcation of the quantum vacuum state may offer an observable signature of Lorentz violation at 
short wavelengths. Indeed, we found that in de Sitter inflation this modiflcation of the vacuum state gives rise to an 
extra overall factor for the super-horizon power spectrum of gravitational- wave background generated from quantum 
fluctuations. 

Since, from a phenomenological point of view, the energy scale of Lorentz violation k Afpj/^/7 can be as low as the 
Hubble parameter during inflation we expect that inflationary cosmological perturbations may open a useful window 
to the physics of Lorentz violation. 



12 
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Appendix A: Weyl-squared actions 



The Weyl-squared action 5w2 — J d'^x \/~g Cabcd Cefgh 1°'^ 1^^ j'^^ u"^ discussed in the main text is one of the 
three invariants consisting of quadratic Weyl tensor fully contracted with ^ab and Ua ; the other two are 



wo 



d^X ^ Cabcd Cefgh 7"' l"^ l'' 7'"' , ^W4 = / d^X Cabcd Cefgh l"' u" u" 



They are such that 



4S^ 



W2 



'WO 



45w4= d^X^CabcdC 



bed 



(Al) 



(A2) 



hence the coefficient —7/2 in (jl.ip and 4-27 in (j2.1l) . 

A useful fact is that these actions are invariant under a conformal transformation, gab = gab ■ the Weyl tensor 
is invariant, C^bcd = C'^bcd , whereas the geometrical quantities associated with the chronon field are transformed 
as Ua = ^ Ua and jab = lab , respectively. Then it can be checked that S'wo , *S'w2 and S'w4 are all conformally 
invariant. 

The expansion of these actions to second order in gauge-invariant perturbations around FLRW spacetimes are (see 
main text for definitions and with = ^ — $ and [3 ) : 



drj d^x 



(A3) 



Appendix B: Derivation of the equations of motion 



To compute the variation of our ghost-free Weyl action with respect to the metric, it is useful to employ an ADM-like 
formalism. To begin with, we express the action as 



(Bl) 



where 



Ua=NdaX, N = 



SW2^ I d'^xV^ Cabcd Cefgh I'^'l'^l^'Wu'^ = j d"" X Wabc W""'' , 

7ab= gab-O-UaUb, a = UaU'', Wabc = lb'' Ic^ u'^ Cdef g ■ (B2) 



1 



VcT daX d'^X ' 



Note that W[ab]c = Wabc , u°- Wabc = Wbca = , and Wab'^ — . We assume a never vanishes so that the lapse 
function N remains finite. We also define the extrinsic curvature of the constant-x surfaces by, V being the covariant 
derivative, 



Kab = la'^lb'^'^dUc 



The Codazzi relation gives 

Wabc ^ la'' lb" Ic^ 



D-2 



V[,(i^el/-^7e]/) 



(B3) 



(B4) 
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where K = 7"'' Kab ■ In order to avoid copious appearances of and Ua , we introduce the following notations: 

7cSa'7/ • • • VdTe../- ^ D,Ta...^ - , u'^ja'lf^ ' • • VrfTe../- ^ V^T^.J - , u«T„...*- ^ T^...''- (B5) 
and so on. With this, we can write 

2 

Kab = DaUb , Wabc = 2£'[a-K'b]c - _ ^ -D[a(-ft'6]c " K 7b]c) • (B6) 

It should be noted that the tensor Wabc only contains one derivative along the direction of Ua ■ 

Let us concentrate on the variation with respect to the metric tensor, gab — >■ Qab + Sgab ■ Define the variables 

A=^u''u^Sgab, Ba = -ia^u''5gbc, Cab=^la''lb'^5gcd- (B7) 

Then we have e.g. 

55(a6 = 2C„b + 2(7U(„B(,)+2crAu„U6, 5^f^ = C + A, (B8) 
where C = 7"^ Cab ■ We also have 

5ua = Aua, 5u" = -Au'' - B'' , 5^ab = 2Cab + 2aU^aBb), 5^a^ = aUaB\ Sj"^ = -2C"\ (B9) 

and 

SKab = V^Cab + 2K^a' Cb)c " iV'^ D^a{N B^)) -KabA + 2a U^a Kb)" B, . (BIO) 
The variation of the Weyl action with respect to the metric can hence be computed, up to divergences, as 

SSw2 = l fd''x^B''>'Sgab, (Bll) 



2 , 
where 

B"^ = 4V„[A^-^ De{N W"'^"''^ )] + 4N-^ Dc{N W^"^^) K - AN'^ D^{N W"^'^") K/^ - AN''^ Dc{N VK«=(»I<*I) iifrf'") 
- ADd[N-^ Dc{N W"^'^")] u^^ - ADd[N-^ D^iN T4^^("l'^l )] i/) + Aa N'^ DdN W"^") K^e v!' 
+ 4Ar-i Dc{N K^^d) + AN-^ Dc{N W"^^" K'>\) + AN'^ D,iN W'^^'^''^ K^a) 

(B12) 

It can be shown that B"'' tensor is traceless and free from the trace part of Kab ■ 

A most important property of i?"** is that it contains derivatives Vu along Ua up to second order only (in the first 
term which is of the type VuDW with W ~ DVu7); hence the eoms remain second order in timelike derivatives 
whenever Ua is timelike. Note also that the other fourth order terms in Bab are all of the type DDW which are 
first order in time and third order in space when Ua is timelike. Because of this distinction between space and time 
derivatives Bab is not Lorentz covariant. 

The variations of the quantities with respect to x given by 

dua = Va , Sjab = -2a U(^a H) whcrc Va = N -jj' dbSx ■ (B13) 
The variation of the Weyl action with respect to x is then computed, up to divergence, as 

6Sw2 = -J d^'x^VaW Sx , (B14) 

where 

Wa = Aa Wabc C''"^" + 2jab W^de C"""" . (B15) 

Note that we may rewrite as Va^F" = -/V"^ Da{N W^) . 
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